The Basics of Water Waves Theory for Analogue Gravity 
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This chapter gives an introduction to the connection between the physics of water 
waves and analogue gravity. Only a basic knowledge of fluid mechanics is assumed 
as a prerequisite. 

"How dare you roil my drink ? Your impudence I shall chastise !" 

"Let not your majesty," the laiTib replies, 

"Decide in haste or passion !" 

"For sure It's difficult to think In what respect or fashion my drinking here could roil your drink, since on the 

stream your majesty now faces I'm lower down, full twenty paces". 

Jean do la Fontaine 
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FIG. 1: The Wolf and the Lamb by Gustave Dore. 



Introduction 

According to Pierre-Gilles de Gennes, "the borders between great empires are often popu- 
lated by the most interesting groups". Indeed, these people often speak several languages and 
are more open-minded due to cultural exchanges. A wonderful analogy exists between the 
propagation of hydrodynamic waves on a fluid flow and the propagation of light in the curved 
space-time of a black hole. It allows us to test astrophysical predictions such as Hawking 
radiation and the effects of high frequency dispersion on it. It provides new insights in Fluid 
Mechanics thanks to the use of tools and concepts borrowed from Quantum Field Theory 
in curved space-time and vice versa. General relativists speak with hydraulicians and this 
chapter is a testimony of their common language and relationships [1, 2]. 

Here, we provide the general background on water waves propagation for analogue gravity: 
we will try to explain how water waves propagate and how a flow current implies the existence 
of an effective space-time; we will insist on the difference between propagation in deep and 
shallow waters on the dispersion relation; a generalized definition of a horizon will be given 
and which, in the particular case of shallow water, reduces to the usual habit of general 
relativists. 

I. A GLIMPSE OF DIMENSIONAL ANALYSIS 

The equations of fiuid mechanics are known since several centuries but they still defy 
modern physics when we try to understand one of its outstanding mysteries like turbulence. 
As they are non-linear and feature several effects such as pressure, gravity and external forces, 
practitioners have been forced to introduce a very useful way to grasp the relevant effects 
when dealing with a peculiar flow. This technique is the so-called dimensional analysis which 
is often the only rescue procedure to disentangle the relative magnitude of several processes 
at play. The reader will be referred to the book by Barenblatt on scaling and dimensional 
analysis for a thorough introduction [3]. Here, we will construct with simple arguments 
the relevant velocities of propagation of water waves depending on the water depth. Three 
regimes of propagation will be uncovered with corresponding dispersion relations. Then, the 
effect of a current will be added and this ingredient is going to be essential in order to have 
an analogue gravity system. 



Waves are characterized by both their wavenumber k and their angular frequency u. Since 
we are deahng with Newton's laws of motion applied to fluids, the second time derivative 
(namely inertia) will translate in a square term in the angular frequency and the dispersion 
relation has the general form w^ = F{k). Obviously, waves with both positive and negative 
angular frequency are thus described by the dispersion relation {cu = ±y^F{k)). Usually, the 
negative root is dismissed since, when the propagation is free, it is a matter of convention to 
focus either on the right or left-propagating waves. Of course, in the presence of a current, 
the system will no longer be symmetric with respect to space reflection {x — )■ —x) and then 
both types of waves (positive and negative) turn to be important... 

A. Shallow waters 

Let us assume that a train of continuous sinusoidal water waves with wavelength A is 
propagating at the surface of a fluid at rest in a given depth h. We make the strong 
hypothesis that the wavelength is much longer than the depth, that is kh << 1 using the 
wavenumber k = 2tt/X. The fluid vertical extension h has the dimension of a length L. Since 
inertia is balanced by gravity, the gravity field is a relevant parameter and its intensity g has 
the dimension of an acceleration L.T^^. We are looking for the typical scaling of the wave 
velocity. The crests of the water waves propagate with the so-called phase velocity and its 
value c has the following dimension L.T~^. We are lead to the obvious scaling law up to a 
constant term : 

Cshaiiow ~ y gh (1) 

Using the definition of the phase velocity c^ = u/k, it is straightforward to infer the 
approximate dispersion relation for water waves propagating in shallow waters : 

io^ ^ ghk^ (2) 

which is similar to the dispersion relation for light waves in empty flat space-time. Analogue 
Gravity will emerge when effective curved space-time is added as we will see... 

Gravity waves in shallow waters are not dispersive since, whatever their wavelength, they 
do propagate with the same velocity. 



B. Deep waters 

Far from the sea shore or for very short gravity waves, the water depth is no more a 
relevant parameter and the only length scale left is the wavelength of the water waves. So, 
if we assume that kh » 1 and recalling that the wavenumber k has the dimension of an 
inverse length L~^, we find the scaling for the phase velocity in deep waters: 

Cdeep ~ y r i^) 

and the dispersion relation: 

cu^ ^ gk (4) 

Newton derived this scaling in his Principia by applying the Galileo formula for the 
period of oscillation of a pendulum T ~ yQg to the water waves. Indeed, if the length of 
the pendulum / is replaced by the wavelength, we do recover the same scaling. 

Gravity waves in deep waters are dispersive since longer waves propagates faster than the 
short ones. 

C. Arbitrary water depth 

Without doing more calculations, we can anticipate that the general dispersion relation 
for waters on arbitrary depth will write according to the following form : 

tu"^ =gk* H{k) (5) 

bearing in mind that the following asymptotic limits must be fulfilled : \m\kh~¥oo H{k) = 1 
and limfc/i^o H{k) = kh. It can been shown that one has rigorously H{k) = tanh (kh) after a 
lengthy calculation that we will avoid to the reader (see the books by Mei [4] or Dingemans 
[5] for a mathematical demonstration). 

D. The capillary length 

Water is made of molecules. In the bulk of the liquid, every molecule is surrounded by the 
same number of neighboring molecules whereas at the boundary with another substance (like 



air at the free surface of sea water) , a water molecule is also surrounded by gas molecules and 
the resulting imbalance in the chemical interactions results in a pressure difference between 
the liquid and the gas. The Young-Laplace law states that the pressure jump is proportional 
both to the local curvature of the interface and to a phenomenological coefficient 7 named 
surface tension which is a property of both media : 

ApY-L = iC (6) 

where C is the curvature of the free surface and has the dimension of an inverse length 
(roughly the radius of curvature). 

When the Young-Laplace pressure Apy-i is balanced by the Pascal pressure drop App = 
pgh due to the static gravity field, a scaling law is deduced easily for the so-called capillary 
length which is the typical size on which surface tension effects are acting : 

ApY-L ~ r = ^Pp ~ P9^^ C^) 



Hence, the capillary length writes: 




lc=^- (8) 

pg 



Let us introduce the effective gravity field g* induced by the capillarity and its dispersive 
scaling law in terms of the wavenumber k: 

S-^^-^'-k^ (9) 

The dispersion relation for water waves taking into account the effect of surface tension 
becomes: 

u^ = {g + g*)ktanh{kh) (10) 



that is [4]: 



u^ = I gk + -k''] tanh (kh) (11) 



II. LONG WATER WAVES ON A CURRENT AS A GRAVITY ANALOGUE 

In a letter to H. Cavendish in 1783, Reverend John Michell introduced the concept of 
what is named, in modern physics, a "black hole" (he was inspired by the corpuscular theory 
of light by I. Newton) [6]: "// the semi-diameter of a sphere of the same density as the Sun 
in the proportion of five hundred to one, and by supposing light to be attracted by the same 
force in proportion to its [mass] with other bodies, all light emitted from such a body would 
be made to return towards it, by its own proper gravity". According to him, this situation 
occurs when the "escape velocity" of a massive particle is equal to the velocity of light. Then, 
Pierre-Simon de Laplace introduced the term "etoz/e sombre" (dark star) in his "Exposition 
du Systeme du Monde" in 1796 to denote such an object. Laplace is also well known for 
having proposed in 1775 an analytical model to describe standing water waves in shallow 
water and for having derived the related dispersion relation [7]. 

Recently, Schiitzhold & Unruh derived the equation of propagation of water waves moving 
on a background flow in the shallow waters limit [8] and this equation describes also the 
behavior of light near the event horizon of a black hole. Indeed, under the impulsion of 
the seminal work by Unruh [9], there has been more and more interest for analogue models 
in general relativity in order to understand the physics of wave propagation on an effective 
curved space-time. Several systems exhibit a so-called "acoustic" metric similar to the metric 
describing a black hole when a wave is moving in a "flowing" medium [1, 2]. 

Here, we will reproduce the derivation of Schiitzhold & Unruh with some details for 
pedagogy [8]. Hence, we consider the propagation of small linear perturbations of a free 
surface between water and air in the presence of an underlying current. The current is 
uniform in depth z, time-independent and varies slowly in the longitudinal direction x. We 
are in the so-called WKBJ approximation such that the wavelength is smaller than the 
typical length on which the current varies (A << -s§^). U = v^ will denote the background 

dx 

flow current whereas v will stand for the velocity associated to the propagation of waves. 
The liquid is inviscid, its density is constant (p = const) and the flow is incompressible. 
We have the following equations of motion : 

- the continuity equation which comes from the incompressibility condition: V.v = 

- the Euler equation : , __ „ 

= v+ v.Vv = ^ + g + - 12 

dt p p 
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with p the pressure, g = —gCz the gravitational acceleration and f = — pV||K" a horizontal 
and irrotational force in the x direction (||) driven by the potential V" which is at the origin 
of the flow. 

Since the flow is assumed to be vorticity-free V x v = (a crucial feature of water waves 
propagation), one has (v.V)v = (V x v) x v + ^V(f^) = ^V(f^) with v = V0 where (j) 
stands for the velocity potential. The vectorial Euler equation reduces to the simpler scalar 
Bernoulh equation 

4> + AV4>f = ---9z-v^^ (13) 

z p 

The boundary conditions are such that: 

- in z = 0, the vertical flow velocity must be null, i.e. ^-'"(2; = 0) = 0. 2; = is by definition 
the bottom depth. 

- the height variations of the fluid are determined by the very same velocity but computed 

on the free surface: ,, 

y^(^z = h) = — = h+{^.V)h (14) 

(JjL 

where ^ is the velocity of a point on the air-water interface. 

- the relative pressure with respect to the atmospheric pressure on the free surface cancels 
by definition : p{z = h) = 0. 

Let us consider a velocity perturbation 6v of the background fiow v^ with a corre- 
sponding vertical displacement 6h. We assume the background fiow v^ to be stationary, 
irrotational and horizontal: V_lVb = 0, v^ = v^ — )> Vy.v^ = 0. 

The Bernoulli equation gives: 

Ivl = -J-9^- y (15) 

where ps follows Pascal's law for static pressure distribution in the water column 
PBiz) = pg{h- z). 

We assume that the velocity perturbation 6v is also curl-less: hence, we can define a 
perturbed velocity potential 6(f). Using again the Bernoulli equation, we get: 



II Sv 

50 + v!j.V||50= — - (16) 

By taking into account the condition pb{z = h) = and using the expression for pb, we 
obtain the boundary condition for the pressure fluctuation at the free surface 6p: 6p{z = 
h) = gp6h. 

The same procedure apphes to the vertical velocity, using the following conditions (14) 
and v^{z = 0) = 0: 6v^{z = 0) = and 6v^{z = he) =6h+ (v|.V||)(5/i • 

We now develop the velocity potential 6(f) using a Taylor series: 

°° z^ 
(50(a;, I/, z) = ^ —6(^(n) {x, y) (17) 

n=0 

The boundary condition v-^{z = 0) = 0, implies (50(i) = 0. With the continuity equation, 

we find : 

V550(o) + (50(2) + ... = O (18) 

We assume that the depth h is much longer than the wavelength A of the free surface 

perturbation. Hence, the higher-order terms in the Taylor expansion are suppressed by 

powers of h/X <^ 1 since we have Vn = (9(1/A^). Keeping only the two lowest terms in the 

Taylor series, we get: 

Sv^{z) =Vx50 

= Vx (50(0) + ^50(2)) (19) 

= Z6(f)(^2) 

At the free surface z = h and using (18), we then find: 

Sv^{z = h) = -hVp(j)(^o) (20) 

In order to find a wave equation for (50(o), let us take the partial derivative with respect 
to time t of the equation (16), then, with the help of the boundary conditions 6p{z = h), 
6v-^{z = h) and of the equation (20), we subsitute 5h with an equivalent expression. Finally, 
we end up with: 

9,2((50) + 2(v|.V||)9t(50) + (vj, ® vjj - gh)V^{6(j)) = (21) 

that is the so-called Beltrami-Laplace equation: 



0(50(0) 



\J-9 



d^^iV^g^'dM^o)) = 



(22) 



provided that we identify g^^^ as an inverse metric expressed in the following matrix form: 



/ 



fiv 



.11 \ 



' B 



1|2 



(23) 



\ A '■■ ^B - ghi J 



Using g^^g^a = 5^ we get the so-called acoustic metric g^y in its typical Painleve- 
Gullstrand form: / 



1 



gh-^B ■ 



\ 



(24) 



/ 



V v^ ; -1 

with c = ^/gh the velocity of water waves in shallow water which is the analogue of the 
velocity of light. 

It is straightforward to show that the dispersion relation associated to equation (21) is: 



2 _ „2u2 



{u - k.vy ^ c'k 



(25) 



which describes the propagation of long water waves on a given flow U = v^ where we 
insist on its approximate nature {kh << 1). The flow induces a Doppler shift of the angular 
frequency: we refer the reader to the hydrodynamics literature where the effect of a current 
on water waves has been discussed extensively [10-15]. One speaks of blue-shifting (red- 
shifting) when the current encounters (follows) the waves and the wavenumber increases 
(decreases). When the flow vanishes, the Beltrami-Laplace operator reduces to the usual 
d'Alembertian operator. 

The dispersion relation is solved graphically in Figure 2. The so-called transplanckian 
problem arises when U = —c that is when the wavenumber of the positive solution (in 
green) diverges at -|-oo, disappears and then reappears as a new diverging negative solution 
(in blue) at — oo for increasing modulus of the flow velocity. 
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kh « 1 



u = o 



U<Q 




uj = n{k) = Uk±ck 




u <0 



f/<0 



FIG. 2: Graphical solutions of the dispersion relation (w — k.U) ~ c k . co = ^{k) is plotted as 
a function of k for increasing modulus of the background flow U < 0. The conserved frequency 
u) is the horizontal red dotted line. The green (blue) color corresponds to the positive (negative) 
branches. 

III. FLUID PARTICLES' TRAJECTORIES 

In this part, we recall (without demonstration) some classical results from water waves 
theory on the trajectories of the fluid particles beneath a water wave [4]. For small wave 
amplitudes {ka << 1), the non-linear terms of the Euler equation and of the boundary 
conditions can be neglected. In 1845, G.B. Airy derived within this approximation the fluid 
particles' trajectories compatible with the following dispersion relation valid for pure gravity 
waves without a background flow for a given depth: 



u 



gk tanh (kh) 



(26) 



We denote z' = the mean water depth of the free surface without wave. Let us consider 
the following perturbation with respect to rest: 



z' = ri{x, t) = asin(a;t — kx) 



(27) 
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Airy computed the resulting velocity profile: 



and 



, , , cosh(k(z + h)) . . , , _^. 

u{x,z,t) = au , ,, , ■ smlut — kxj (28) 

sinh{kh) 



. I , sinh(A;(z + /i)) , , -. ,^^. 

w[x,z,t) = au , ,, , , cosiut — kxj (29) 

simi[kh) 



where u and w correspond to the projections of the perturbation velocity in the horizontal 
and vertical directions. 

Under the hypothesis of small displacements, one deduces the horizontal motion of fluid 

particles: 

^^/ / N ^r cosh(A;(2;' + h)) , , , _„, 

X a;,/,t -Xo = -a w, ,x ^^ cos a;t - fca: (30 

smh(A;/i) 

as well as the vertical motion: 

ry, I \ ry Sinh(A;(z' + /l)) . / , X .„,s 

Z(x, z , t) — Zq = a , ,, , , sm(wt — kx) (31) 

smh(A;/i) 

In deep waters (far from the sea shore for example), the fluid particles' trajectories are 
circular with radius R. In shallow waters (close to the beach!), the trajectories flatten and 
the particles follow an ellipse of semi-axis A and B for respectively the horizontal and vertical 
motions. In practice, practitioners distinguish three zones : 

• the deep water case ihl\ > 1/2) : 



A ~ 5 ~ ae"' (32) 

The trajectories are circles of radius R r^ A r^ B which decreases exponentially with 
the depth z' . 

the intermediate case (1/20 < h/\ < 1/2) : 



and 



A _ ^ cos\i{k{z' + h)) 



sin\i{k{z' + h)) 
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[a] Deep waters 
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(&) Shallow waters 



FIG. 3: The flow generated beneath a surface waves: (a) deep water case; (b) shahow water case. 

The trajectories are eUipses whose semi-axes diminish with depth. The decrease is 
slower than the exponential one of deep waters. 

• the shallow water case {h/X < 1/20) : 



A 



kh 



and 



B 



ak{z' + h) 
kh 



(35) 



(36) 



The trajectories are ellipses whose major semi- axis A is independent of the water depth 
z' and whose minor semi-axis B decreases linearly with z'. On the bottom {z' = —h), 
B cancels and the trajectories become a horizontal oscillation of amplitude A. 
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In the presence of a current U, the previous expressions for the velocity field keep the 
same form provided the dependence of the amplitude (and not of the phase) with the angular 
frequency u is replaced by the relative angular frequency oj' = u — k.U [10-16]. Of course, 
u becomes u' + U whereas w' is invariant. The particles' trajectories are thus similar to 
cycloids whose amplitude decreases with the water depth for waves following the current 
[16]. 

IV. A PLETHORA OF DISPERSIVE EFFECTS 

One of the salient effects of analogue gravity is the possibility to solve the transplanckian 
problem thanks to the introduction of dispersion close to the horizon of an artificial black 
hole. As a matter of fact, a major drawback of the original calculation by Stephen Hawking 
of the black hole radiation is the necessity for the field to have a wavelength which goes to 
zero as one gets close to the event horizon. Water waves provide several regularization scales 
in a cascade such as the water depth, the capillary length or even a viscous scale in order 
to cope with a diverging wavenumber by counter-acting the continuous blue-shifting of the 
flow... 

Let us consider the propagation of gravity waves (without surface tension 7 = for the 
moment) on a linear shear flow U{z) = Uq + Viz with constant plug flow Uq and constant 
vorticity Vt. Here, one assumes that both the bottom depth and the flow velocity vary slowly 
such that ^ >> A and ^ >> A. The dispersion relation between the frequency ^ and 

dx dx 

the wavenumber k writes either with its implicit expression due to Thompson [17]: 

(a; - Uokf = [gk - n{uj - kUo)] tanh (kh) (37) 

or with its explicit expression due to Biesel [18]: 



n 




u = Uok — — tanh (kh) ± W I — tanh (kh) | + gk tanh (kh) (38) 

With surface tension, Huang has derived recently the following dispersion relation with 
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its implicit expression [19]: 

(w - kUof = [gk + -P - n{uj - kUo)] tanh (kh) (39) 

which can be written exphcitly according to Choi [20] in the form: 



uj = Uok tanh (A;/i) ± W (- tanh (A;/i) j + ( ^A; + -A;3 j tanh (A;/i) (40) 

It is interesting to notice that the "relativistic" dispersion relation {u — U.kY = c^k"^ 
is recovered in the long wavelength limit kh '^ 1 whatever is the dispersive correction. It 
obvious when dealing with the surface tension since the capillary length is smaller that the 
long wavelength. It is less obvious for the dispersive effect of vorticity. Indeed, the long 
wavelength approximation of the Biesel's dispersion relation writes: 



oo^Uok-^k± Jghk^ + ^P (41) 

Fortunately, it can be transformed into the usual dispersion relation associated to the 
acoustic metric {u — U' .k)"^ = d'^k'^ provided one introduces renormalized flow and waves 
velocities U' = Uq - VLh/2 and c' = ^gh + Vf^K^/A. 

Assuming a uniform flow in the vertical direction (i7 = 0), the dispersion relation becomes 
[10-15]: 

(w - Ukf ~ (gk + -kA tanh(A;/i) (42) 

• In the shallow water limit kh « 1, 

{uj - Ukf ~ ghk^ +(—- —\ k^ + 0(A;6) (43) 

the dispersion relation is identical to a BEC-type phonons spectrum: 

{u - Ukf ~ c^k^ ± c^^k* (44) 
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with the corresponding "heahng length": 



e 



' 3 



(45) 



An interesting observation is that the superluminal correction can have a negative sign 
in contrast to the BEC case if the capillary length is less than /i/v3 or even null... 

• In the deep water limit kh » 1, the dispersion relation looses its "relativis- 
tic/acoustic" branch: 

(w - Ukf ^gk+ -k^ (46) 

Viscosity has both a dissipative (imaginary term) and a dispersive (real term) contri- 
butions to the dispersion relation. By dimensional analysis, it is obvious that the typical 
viscous scale would be of the order of 5 ~ a/^ otherwise known as the Stokes viscous length 
which is the scale of the viscous boundary layer [4] . 

V. HYDRODYNAMIC HORIZONS 

In this part, we propose a generalized definition of a horizon with respect to the usual 
custom in General Relativity. Condensed matter horizons and here, hydrodynamic horizons 
lead to a dispersive-like definition. What is a Horizon ? The word horizon derives from 
the Greek ^^opL^uu nvuXoq" {horizon kyklos), "separating circle", from the verb "opt^u" 
(horizo), "to divide, to separate", from the word "opoq" (oros), "boundary, landmark". In 
the Fable of Jean de la Fontaine recalled at the beginning of this chapter, will the Lamb 
be right to argue against the Wolf that the waves he creates as he is drinking at the river 
border will not climb against the current and reach the Wolf ? Will a frontier separate the 
Lamb from the Wolf: will a horizon form ? Will the position of the frontier depend on the 
period of the waves: will the horizon be dispersive or not ? 

A. Non-Dispersive Horizons 

The analogy between the propagation of light in a curved space-time and the propagation 
of long gravity waves on a current features the so-called "acoustic/relativistic" dispersion 
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relation {u — k.U) ~ c^A;^ as a common characteristic for both systems assuming kh << 1 
and without surface tension. A simple dimensional analysis of it: 

u^ ^ U^e ^ ghk^ (47) 

allows to infer scaling laws for the wavenumber: 

k ^ -^ (48) 

Vgh 



and the blocking velocity: 

U ^ ^^ (49) 

which we confirm by solving the "relativistic" dispersion relation as a polynomial in k: 

U + y'gh 



implying the transplanckian problem {kh — )■ oo) when the blocking velocity [U* in modulus) 
of long gravity waves matches the current flow: 



Uh = -Vah (51) 

This non-dispersive definition of a horizon in hydrodynamics corresponds to the definition 
of General Relativity where the pure temporal matrix element of the Painleve-Gullstrand 
metric is cancelled: 

^00 = (52) 

leading to: 

U = -c=-^^ (53) 

Here, it is crucial to understand that the relativistic horizon hides in fact three intricate 
horizons (Figure 4: blocking velocity U* versus the wave period T): a group velocity horizon 
(cg = ll = f/ + c = 0), a phase velocity horizon (c,^ = | = f/ + c = 0) and a negative 
horizon (or negative energy mode horizon): negative relative frequencies u — k.U < 
can appear. This last fact implies that Stimulated Hawking Radiation can be observed in 
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kh « 1 
>T 



Nq Blocking 



U* = -c 



Blocking 

=> Non-dispersive "horizon" : 
group velocity horizon=phase velocity horizon=negative horizon 

FIG. 4: Phase-Space -U versus T- of the dispersion relation (cj — k.U) ~ c^k"^. 

Classical Physics using water waves and this is one of the major interests of the analogue 
gravity program for the Fluid Mechanics community [21-25]. As soon as there is a phase 
velocity horizon, this one is identical with a negative horizon. Because of dispersion, a phase 
velocity horizon can be absent whereas a negative group velocity horizon can be present (see 
below). 

These "negative energy waves" are well known in Hydrodynamics. Werner Heisenberg 
discovered them in his PhD Thesis on the stability of the plane Couette flow. He showed 
that viscosity can have a destabilizing effect if negative energy waves (also named Tollmien- 
Schlichting waves) are present (at the so-called critical layer corresponding to a phase velocity 
horizon) in a unidirectional non-inflectional plane flow which is normally stable if inviscid 
according to the classical Rayleigh criterion [7, 14]! 



B. Dispersive Horizons 

How is the deflnition of a horizon modifled in the presence of dispersion ? Wave blocking 
is a process where a flow separates a free surface into a flat and a deformed surface. The 
boundary deflnes a "horizon". A wave phenomenon implies the existence of a dispersion 
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relation u = Q{k). At the boundary, the energy flow of the system "waves+current" cancels: 



wave+current n 

^group ~ Q^ ~ ^ 



(54) 



This last criterion will define a hydrodynamic horizon as a group velocity horizon (or turning 
point using WKBJ terminology). Of course, one recovers the non-dispersive definition U = 
— c for an "acoustic/relativistic" dispersion relation. We treat here the simple case of a white 
hole horizon which is the time reverse of a black hole horizon. As previously, the dispersion 
relation for water waves in arbitrary depth is solved graphically (Figure 5). An extremum 
of the function to = il(k) = Uk ± ^^/gkta,nh(kh) corresponds to a horizon. 



c/ = o 



U<0 




FIG. 5: Graphical solutions of the dispersion relation (a; — k.U) = gktanh{kh). The conserved 
frequency uj is the horizontal red dotted line. The green (blue) color corresponds to the positive 
(negative) branches. 

Some scaling laws can be derived in the high dispersive regime where kh » 1: 
• without surface tension. Dimensional analysis leads to: 



uj^ ^ U'^k'^ ^ gk 



(55) 
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that is : 



u 



(56) 



and: 



K UJ 



(57) 



The rigorous mathematical treatment gives [22]: 



Au' 



(58) 



and 



U„ = - 



9 gT 



(59) 



4:00 Sir 

The blocking velocity U* depends now on the incoming period of the water waves 
(Figure 6). Depending on the period, we have either U* = Uh for long waves or 
U* = Ug for short waves. 




T(s) 



LO ^ n(fc) —Uk± y/gk tsinh{kh) 



FIG. 6: Superposed Phase-Spaces -U versus T- for the shallow and deep water cases without surface 
tension. 



with surface tension. Dimensional analysis leads to: 



P 



(60) 



that is: 



and: 
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A:^(^)''' (61) 



U^ {^] (62) 



The rigorous mathematical treatment gives [23]: 

*.-(?)'" (63) 

and: 

U, = -V2 f'^] (64) 

A new horizon (in fact two) appears. Blue-shifted waves and negative energy waves 
can be reflected at a blue horizon and a negative horizon whose common asymptotic 
value is U^ (see [23] for the details and the corresponding chapter in this book). Two 
maxima and a minimum appear in the graphical analysis of the dispersion relation 
(Figure 7). A cusp where the white and blue horizons merge appears in the Phase- 
Space (Figure 8). 

When the water depth changes, the dispersion relation is either {u — f/A;)^ ~ c^k"^ ± 
^2^2^4 £qj, kh << 1 allowing dispersive corrections (only a negative horizon remains 
with the positive quartic correction) or (cu — k.U) = igk + -k^ ) ta,nh{kh) and three 
horizons are observed (Figure 9). h* = 2lc determines the transition depth between 
both behaviors. 
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FIG. 7: Graphical solutions of the dispersion relation {u — k.U) = igk + -k^j tanh(A;/i). The 
conserved frequency uj is the horizontal red dotted line. The green (blue) color corresponds to the 
positive (negative) branches. 



with vorticity. Similar arguments would lead to a new horizon replacing U^ when 
including vorticity Q with qualitatively the same behavior in the limit kh » 1: 



and 



kn 



Uc, 



7 



1/3 



P 



(65) 



(66) 



C. Natural and Artificial Horizons 



In this part, we give some examples of water wave horizons. In hydrodynamics, white 
holes are more usual than black holes whose canonical example is the draining flow in the 
bathtub. A river mouth dying in the sea is a nice case of a natural white hole: the sea waves 
are blocked by the river flow. Figure 10 is an example found by the author when he used to 
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'^ = a(k) = Uk± Jfgk + ^kA 




IUgl = gT/8n 
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FIG. 8: Phase-Space -U versus T- for the deep water case including surface tension. 
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FIG. 9: Phase-Spaces -U versus T- for a changing water depth including surface tension. 
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walk on the Promenade des Anglais in Nice (France). This white hole inspired the following 
studies [21, 23, 24]. 




FIG. 10: A natural white hole of the French Riviera. 



A more controlled white hole in the laboratory was suggested a few years ago by Volovik 
[26]: the circular jump in the kitchen sink (Figure 11). We studied its related Mach cone 
and its dispersive properties in [25]. 

Recently, the author became aware of a biological-induced white hole with interesting 
dispersive properties, namely the whale fluke-print. As a whale swims or dives, it releases 
a vortex ring behind its fluke at each oscillation. The flow induced on the free surface is 
directed radially and forms a oval patch that gravity waves cannot enter whereas capillary 
waves are seen on its boundary (Figure 12). 

Dispersive and non- dispersive horizons are all encountered in nature and can be simulated 
in the laboratory. Dispersion has another intriguing consequence, namely the appearance 
of current-induced zero-frequency waves which appear spontaneously. These so-called "zero 
modes" have no counterpart in General Relativity so far... 
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FIG. 11: A laboratory white hole in the kitchen sink. 




FIG. 12: A "biological" dispersive white hole. 
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D. Zero Modes 

A flat interface can be considered as a wave witli zero frequency u = and zero wavenum- 
ber k = 0. When a spatially varying current is flowing under such an interface, the infinite 
wavelength of the interface can be reflected at a blocking line (creation of a group velocity 
horizon). This process produces a static (w = 0) jump through the interferences between 
the incident wave (flat surface with infinite wavelength k = 0) and the reflected one. This 
explains the formation of the circular [25, 26] and hydraulic [27] jumps. An undulation is 
observed which has a zero phase velocity but a non-zero group velocity and thus withdraws 
energy from the horizon towards infinity. When surface tension is present, the "gravity" 
jump is decorated by static capillary ripples inside the circular jump [28]. 

1. The zero mode (static undulation) for gravity waves 
Two opposite wavenumbers are solutions of the dispersion relation for a zero frequency : 

{O-Ukf = gkt8.nh{kh) (67) 

Clearly, there is no threshold since there is always a solution whatever the velocity of the 
flow. Whatever the water depth, the slightest current flow induces a free surface deformation. 
Let us take the extreme shallow waters limit {kh << 1): 

U^"^ ~ ghk"^ (68) 

The threshold velocity for the zero mode appearance would correspond to: 

'-^ — V yi^ Cphase ^group V / 

that is: 

U , , 

Fr = —= = 1 (70) 

vgh 



in terms of the dimensionless Froude number Fr. This latter constraint is well known in 
Hydraulics as the condition of appearance of the hydraulic jump when water flows over a 
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bump: 

Froude = Fr = —= = = = M = Mach (71) 

y (jH Cp/ifl-se (^ group 

which is similar to the supersonic-subsonic transition of air flows in aerodynamics described 
by the so-called Mach number M [25]. 

Then, one distinguishes in Hydraulics the following regimes: 

• Fr < 1: (a) subcritical-to-subcritical flow over a bump. A group velocity horizon can 
appear but no phase velocity horizon [21]. No hydraulic jump is created but a static 
undulation is observed. The water depth decreases on average over the bump. 

• Fr > 1: (b) supercritical-to-supercritical flow over a bump. The water depth increases 
on average over the bump. 

• Fr = 1: (c) subcritical-to-supercritical flow over a bump. The group and phase veloc- 
ity horizons are the same. A hydraulic jump appears as part of the static undulation 

[24, 27]. 

2. The zero mode (static undulation) for capillo- gravity waves 
If one takes into account the effect of surface tension: 

(0 - Ukf = (gk + -kA tanh{kh) (72) 

One is lead to the existence of a velocity threshold which corresponds to the minimum 

/ \ 1/4 

of the phase velocity with the wavenumber Uj = —\/2 ( — j [23]. 
For the case of the circular jump assuming kh << 1, 

(0 - Ukf = c'k^ + (ll - J^ k* (73) 

the following condition: 

h < ^/Slc ^U >c= ^^ (74) 

implies the existence of static capillary undulations in the supersonic region of the circular 
jump [28]. 
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We have seen how the dispersion relation explains the appearance of a horizon as well 
as the evolution of the wavelength of the converted modes. How does the amplitude of the 
modes evolve ? 

VI. THE "NORM" 

// someone tells you that he knows what E ~ hw means, 
tell him that he is a liar. 

Albert Einstein 
In this final part, we will show how the so-called "norm" used by relativists in order to 
derive the Hawking spectrum is nothing else than the wave action, a pure classical concept. 
In 1905, Albert Einstein pointed out that the four-momentum and the four-wave vector 
transform similarly under a Lorentz boost. This simple remark was fundamental in order 
to infer the existence of the light quantum whose energy is proportional to the frequency. 
The factor of proportionality was the Planck constant and Physicists soon realized that the 
latter constant of nature was measured in units of action. Quantum Mechanics then will 
shortly take its roots in Analytical Mechanics. In the famous 1911 Solvay conference [29], 
Lorentz wondered about the paradoxical behavior of a harmonic oscillator like a pendulum 
whose frequency was made to change slowly with time by reducing its length. Indeed, the 
corresponding quantum behavior of the oscillator would forbid a change in the quantum 
number describing the state of the oscillator since the frequency variation would not be high 
enough to allow transition to another state. Einstein pointed out that both the energy and 
the frequency of the pendulum would change with time but not their ratio as discovered 
by Rayleigh in 1902 [30]. Ehrenfest showed that the ratio of energy to frequency namely 
the action was an "adiabatic invariant" . Adiabatic invariants of a given dynamical system 
are approximate constants of motion which are approximately preserved during a process 
where the parameters of the system change slowly on a time scale, which is supposed to 
be much larger than any typical dynamical time scale. They are the quantities to quantize 
when switching from Analytical Mechanics to Quantum Mechanics. A similar relation was 
discovered later by De Broglie between the momentum and the wavenumber. It should be 
borne in mind that the ratio of the action to the Planck constant is the number of photons 
which is another way to interpret the norm (as we will see) as the number of photons/phonons 
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times the quantum of action in a quantum context. 

Water waves are an example of a classical field. Thus, we can anticipate that the fluid 
system will have a corresponding wave action density defined as the ratio between a mean 
energy density (computed by averaging the instantaneous energy density on a spatial wave 
period) and the wave frequency. This action J = E/u is assumed to be an adiabatic invariant 
(see [31] for a demonstration based on classical field theory). Then, if prime denotes a moving 
frame of reference with velocity v, we must have J' = J for a Galilean boost (recall that 
Einstein dealt with Lorentz transformations applied to light) that is [32, 33] : 

E' E 

- = - 75 

which is valid if and only if we have the following transformations : 

E' = E- v.P, P' = P (76) 

and 



UJ' 



' = 00- v.k, k' = k (77) 



The latter formulae are just the usual Doppler effect whereas the former correspond to the 
change of energy/momentum for a classical wave and NOT a particle [34]. These would 
apply to a quasi-particle that is a collective excitation : phonon in acoustics or ripplons for 
water waves. In the following, momentum and energy would refer to quasi-momentum and 
quasi-energy if not specified. 

It is now obvious that the energy in the moving frame will be given by E' = E{1 — v/c^) 
where c^ = oo/k is the phase velocity in the rest frame [32, 33]. In order to have negative 
energy waves {E' < 0), the so-called Landau criterion must be fulfilled oo — U.k < since 
the energy in the rest frame i5 > is always positive. It is well known that superfluidity 
is lost when negative energy waves are created at the minimum of the roton spectrum [26] . 
It is similar to waves creation (Cerenkov-like effect) by an object in a flowing current U 
perforating the interface between water and air (capillary waves in the front and gravity 
waves in the rear) when min(c0) < U : the phase velocity c^ = u/k features a minimum 
under which no waves are created (cf. Thomson and Helmholtz flshing line as described by 
Darrigol [7] and the corresponding chapter in this book). Let us recall that, in the direct 
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space, the mean energy density (or pseudo-energy) for water waves (without any current) is 
proportional to the square of the amphtude E = l/2pga'^ where a is the amphtude of the 
wave [4]. 

Here, we must be careful when we want to evaluate the wave energy in the moving frame 
of the current because the velocity of the flow is changing with space. Hence, every part 
of the water waves wavelength will be " desynchronized" by the spatial-dependent Doppler 
effect due to the current. That is why Weinfurtner et al. [24] introduced a time shift 
tc = f jjf^ which is reminiscent of Carroll kinematics for classical waves t' = t — vx/c^ and 
x' = X or t' = t — x/Vq where Vq = (? jv is the dual velocity associated to the wavefront 
[34]. The Carrollian time shift writes di! = dt — dx/Vo in differential form where Vq is now a 
function of space x in the water waves problem. The tc coordinate has dimension of time and 
its associated "wavenumber" fc has units of a frequency. The usual convective derivative 
operator dt + U{x)dx becomes dt + dt^ and in Fourier transform space / -|- f^- Then, when 
analysing data in the Carrollian coordinate system, the amplitude of the wave 77 is a function 
of both the normal time t and the Carrollian time tc that is in the Fourier transform space 
vifi fc)- The wave action density in the Fourier space Jwave is by definition the integral on 
the different Carrollian times of the ratio between the Fourier transform of the wave energy 
density (E) and the Fourier transform of the relative angular frequency (w'): 

'Jwave / 7 , 7 ^/c V / 

J J ' Jc 

The expression of Jwave is similar to the Zeldovich formula for the number of photons A^ 
when dealing with plane electromagnetic waves that are not monochromatic {u = ±CL|k|): 

^ = ii/ftEMP±MM)£ (79) 

hSn J ±CL\k\ ^ ' 

where cx, is the light velocity and E(k, t), B(k, t) are the Fourier transforms of the electric 
and magnetic fields [35]. The sign in the denominator comes from the dispersion relation 
of light which features both positive and negative branches. The Zeldovich formula and the 
"norm" used in [24] writes as the ratio between a wave energy (which scales with the square 
of an amplitude) and the wave frequency. The case of acoustics is discussed in [36] following 
the treatment by Landau and Lifschitz [37]. 
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The equivalence between the norm and the wave action density can be formally proven as 
follows. First, wave packets on the free surface of water obey the Beltrami-Laplace equation 
in the Painleve-Gullstrand metric as shown by Schutzhold and Unruh in 2002 (see [8] and 
the corresponding chapter in this book): 

dt{dt<P + f/9,0) + d^iUdtcj) + U^d^cj)) - c^dl<P = (80) 

where is the velocity potential fluctuation. The complete velocity potential featuring both 
the waves and the background flow U is such that its space derivative is by definition the 
flow velocity. 

We can expect the conservation of two quantities due to the invariance of the correspond- 
ing action under 1) the transformation — )■ e*"0, a constant, and 2) time translation (for 
time- independent U). The former invariance gives conservation of the Klein-Gordon norm 
(as demonstrated elsewhere in this book): 

iV = 2^ / dx [0*(9i0 + Ud^<P) - 0(9^0* + Ud^*)] , (81) 

whereas the latter gives conservation of (pseudo-) energy. For wave packets confined to a 
region where the fiow velocity U is constant, the norm (81) can be written in /c-space in 
terms of the Fourier transform 0(A;) as: 

N=- dk{u-Uk)\(t){k)\^, (82) 

'^ J -oo 

The typical interpretation of the Zeldovich formula is that it is a positive quantity: the 
number of photons. However, the Klein-Gordon norm used but the relativists is either 
positive or negative. Then, the Zeldovich formula encodes in general a different information 
than the Klein-Gordon norm. The latter counts the amount of charge, that is why for real 
fields it is zero. The complex solutions however do have charge. In fact, strictly speaking, 
the last equation is not correct as for a single k there can be modes with ±|A;|. Thus, apart 
from the integral there should be a sum in positive/negative branches: 

N = ^^ j dk{c\k\\au\'' - c\k\\h\''). (83) 
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with a^, bk the corresponding Fourier coefficients. 

The Zeldovich formula (without the negative sign in the denominator) does not give zero 
for real fields. It really corresponds to^ J (iA;(c|A;||afcp). It does not contain the second term 
that in the case of real fields (ofc = bk), would combine to yield a total zero (this is the case 
also for photons). The point is that what corresponds to the norm is not one of the pieces 
individually but the addition of the two. 

In general, the norm scales like the integral over the wavenumber of the amplitude square 
of the Fourier transformed velocity potential times the relative frequency in the moving 
frame. Hence, the norm scaling is A^ ~ J dk{u — Uk)(f)'^ [24]. However, because the velocity 
potential is related to the free surface deformation rj by the Bernoulli equation d(f)/dt+gri = 
(here, without a flow to simplify), it follows that the velocity potential scales like ~ gfj/{u— 
Uk) in the Fourier space [8]. We conclude that the norm behaves like N ^ J dkg'^ff/{uj — Uk) 
as the wave action that is as the ratio between the square of the amplitude (the energy) and 
the relative frequency. 

The norm is strictly conserved. Is this the case for the wave action ? In the fluid mechanics 
literature, the wave action is the solution of a conservation equation which replaces obviously 
the conservation of energy of a closed system. Here the system is open since the waves 
interact with the flow and do exchange energy. Bretherton and Garrett have shown that 
the wave action conservation writes in the so-called WKBJ regime where the flow velocity 
varies on a length scale much larger that the wavelength [38, 39]: 

where u' = oo — Uk and Cg is the total group velocity including the background flow. Accord- 
ing to Bretherton and Garrett, ^^because E' is an energy density, it is not constant down a 
ray, even if wave energy is conserved. However, in a time dependent and/or non-uniformly 
moving medium, u' varies along a ray. If E' /u' is the wave action density, total wave action 
is conserved, whereas total wave energy is not" . For a stationary process, we deduce that 
rj^Cg/{oo — kU) = const since the energy density in the moving E' is proportional to the square 
of the interface deformation rj (as in the rest frame without current). As a consequence, the 
amplitude diverges to infinity if one gets close to a turning point where the group velocity 
vanishes and where the WKBJ approximation is no longer valid. Dispersion enters the game 
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to avoid such a caustic. 

The change of wave action J = E/u of a slowly modulated oscillator is exponentially 
small in the non-adiabatic parameter {u / (rate of change of the medium properties)) : a 
mathematical theorem due to Meyer in 1973 [40]. For the linear pendulum of Rayleigh with 
a varying length, the rate of change is directly the inverse of the time lapse. Here, with water 
waves on a non-uniform flow, the property is the velocity U and its typical rate of change 
is its space gradient (the so-called surface gravity in General Relativity) whose dimension 
is the one of a frequency: Jacobson and Parentani defined the surface gravity as a local 
expansion rate seen by a freely falling observer when he crosses the horizon [41]. One is 
tempted to extrapolate the following behavior for the change of wave action as the waves 
propagate against the flow: 

A J = Jo exp ( —constant— j (85) 

where k = dU/dx is the surface gravity for water waves. 

The change of wave action would be very similar to the famous Hawking spectrum [2] : 

^ = exp (-2n-) (86) 

The fact that the bogoliubov coefficient behaves as an exponential has been discussed by 
Jacobson [42]. 

Let us introduce the following dimensionless numbers u / (rate of change of the medium 
properties) with names of distinguished physicists : 

— = Hawking number = J^ (87) 

K 

Unruh number = ^„ (88) 



dx 



The validity of the WKBJ inequality ^ >> A can be reassessed thanks to the Unruh 

dx 

number close to the horizon. As a matter of fact, U* ~ g/u and A* ~ g/oo'^ then ^„ ^ 0(1): 
close to the horizon, the WKBJ approximation breaks down. If ^„ >> 1, then the process 
is adiabatic. In order not to have a vanishing spectrum, ^„ ^ 0(1), then vacuum radiation 
a la Hawking-Unruh is a non-adiabatic process [43]. The case ^„ << 1 would imply a too 
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small frequency, hence the amplitude of the energy spectrum (which scales with the cube of 
the frequency in 3D) would vanish. 

Conclusion 

This rapid tour of the field of analogue gravity through the prism of water waves theory 
has broadened our definition of a horizon and has deepened our understanding of the concept 
of norm as used by relativists. In a related chapter of this book, we study experimentally 
the influence of surface tension and the associated dispersive horizons. Moreover, we try to 
answer to the question "what is a particle close to a horizon? " . 
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